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Abstract. We present a method of localised control of chaos in Hamiltonian systems. 
The aim is to modify the perturbation locally by a small control term which makes the 
controlled Hamiltonian more regular. We provide an explicit expression for the control 
term which is able to recreate invariant (KAM) tori without modifying other parts of 
phase space. We apply this method of localised control to a forced pendulum model, 
to the delta-kicked rotor (standard map) and to a non-twist Hamiltonian. 
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1. Introduction 

Controlling chaotic transport is a key challenge in many branches of physics like for 
instance, in particle accelerators, free electron lasers or in magnetically confined fusion 
plasmas. One way to control transport would be that of reducing or suppressing chaos. 
There exist numerous attempts to control chaos (see Refs. [TJ Ej for a rather extended 
list of references). Most of the methods for controlling chaotic systems is done by 
tilting targeted trajectories. However, for many body experiments like the magnetic 
confinement of a plasma or the control of turbulent flows, such methods are hopeless 
due to the high number of trajectories to deal with simultaneously. For these systems, 
it is desirable to control transport properties without significantly altering the original 
system under investigation nor its overall chaotic structure. Here we focus on another 
strategy which is based on building barriers by adding a small apt perturbation which 
is localised in phase space, hence confining all the trajectories. 

The main motivations for a localised control are the following ones : Very often the 
control of a physical system can only be performed in some specific regions of phase 
space. This is in particular the case in thermonuclear fusion devices where the electric 
potential can only be modified near the border of the plasma. For some purposes it is 
sometimes desirable to stabilize only a given region of phase space without modifying the 
major part of phase space in order to preserve some specific features of the system. This 
method can be used to bound the motion of particles without changing the perturbation 
inside (and outside) the barrier. Also, using a localised control means that one needs 
to inject much fewer energy than a global control in order to create isolated barriers of 
transport. 

In this article, we consider the class of Hamiltonian systems that can be written in 
the form H = H + sV i.e. an integrable Hamiltonian H (with action-angle variables) 
plus a small perturbation eV. The idea is to slightly and locally modify the perturbation 
and create regular structures (like invariant tori) : The aim is to devise a control term 
/ such that the dynamics of the controlled Hamiltonian H c = Hq + eV + / has more 
regular trajectories or less diffusion than the uncontrolled one. For practical purposes, 
the control term should be small with respect to the perturbation eV, and localised in 
phase space (i.e. the subset of phase space where / is non-zero is finite and small). 

In Refs. [3J IU , an explicit method of control was provided in order to construct a 
control term / of order e 2 such that the controlled Hamiltonian Hq+eV + / is integrable. 
The main drawback of this approach is that the control term has to be applied on all the 
phase space. Here we provide a method to construct control terms / of order e 2 with a 
finite support in phase space, such that the controlled Hamiltonian H c = H + eV + f 
has isolated invariant tori. For Hamiltonian systems with two degrees of freedom, these 
invariant tori act as barriers in phase space. For higher dimensional systems KAM tori 
act as effective barriers of diffusion. 

The main result of the paper is stated as follows : For a Hamiltonian system 
written in action-angle variables with L degrees of freedom, the perturbed Hamiltonian 
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is H(A, 0) = uj ■ A + V(A, 9) where (A, 6) e M L x T L and a; is a non-resonant vector 
of M. L . We consider a region near A = and the perturbation V has constant and 
linear parts in actions of order e, i.e. V(A, 0) = sv(0) + ew(6) ■ A + Q(A,6) where 
Q is of order 0(||A|| 2 ). We notice that for e = 0, the Hamiltonian H has an invariant 
torus with frequency vector uj at A = for any Q not necessarily small. The controlled 
Hamiltonian we construct is 



where fHs a smooth characteristic function of a region around a targeted invariant torus 
(the size of its support is of order e). It is sufficient to have f2(||A||) = 1 for ||A|| < e. 
For instance, Q — 1 would be a possible and simpler candidate, however representing 
a long-range control. We notice that the control term / we construct only depends on 
the angle variables and is given by 



where T is a linear operator defined below as a pseudo- inverse of uj ■ de. Note that / 
is of order e 2 . For a sufficiently small perturbation, Hamiltonian has an invariant 
torus with frequency vector close to uj. After proving this result, we check numerically 
that the controlled Hamiltonian is more regular than the uncontrolled one, i.e. the 
invariant tori of the controlled Hamiltonian persist to higher values of the amplitude of 
the perturbation than in the uncontrolled case. 

In Sec. we explain the theory of the localised control of Hamiltonian systems 
and in particular we prove Eqs. JIJ)-(J2J)- In Sec. El we give some applications of the 
localised control on the following models: a forced pendulum Hamiltonian, the delta- 
kicked rotor (standard map) and a non-twist Hamiltonian model. For these systems, 
we show numerically that the localised control is able to create isolated invariant tori 
beyond the values of the parameters for which there are no invariant tori in the absence 
of control. 

2. Localised control for Hamiltonian systems 

2.1. Global control : toward a localised control theory 

We first recall the global control theory as explained in Refs. E] in order to define 

the main operators that will be used for the localised control. 

Let us fix a Hamiltonian Hq. We define the linear operator {Hq} by 



H C (A, 6) = uj ■ A + V(A, 6) + /(0)O(||A||), 



(1) 



f(0) = V{0,0)-V(-Td e V{p,0),0) 



(2) 



{H }H = {H ,H}, 



where {•, •} is the Poisson bracket. The operator {-f^o} is n °t invertible, e.g., {H }H Q = 
0. We consider a pseudo-inverse of {-f^o}; denoted by T, satisfying 

{#o} 2 r = {H }. (3) 
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If the operator T exists, it is not unique in general. We define the resonant operator 7Z 
as 

K = l- {H }T, (4) 

We notice that Eq. (J3J) becomes {H }7Z = 0. A consequence is that any element 1ZV is 
constant under the flow of Hq. 

Notation : In what follows, we will use the notation ab for an operation between 
a and b which can be vectors or covectors. For instance, if a is a covector and b a 
vector, ab is the usual scalar product a • b. If a is a vector and b a covector, ab is 
the matrix whose elements are [ab]^ = a$ . For a vector a and a matrix M, Ma is a 
vector. In the same way, if a is a covector, aM is a covector. Also we denote a' the 
matrix c^a with elements [cfoa]^ = da? jdQ % . For clarity we also denote d the operator <9#. 

Let us now assume that Hq is integrable with action-angle variables (A, 0) e 
M. L x T L where T L is the L-dimensional torus. Here, A is a L-dimensional vector 
and 6 is a L-dimensional covector. The Poisson bracket between two functions H and 
V is given in the usual form 

_ OHdV _ dVdH 
* ' ^ ~ 8A~d0 ~ 8A~d0' 

We assume that Hq is linear in the actions variables, so that Hq = u>A, where the 
frequency vector uj is any co- vector of M. L . In this paper, we assume that u> is non- 
resonant, i.e. there is no vector k e 1 L \ {0} such that u>k = 0. The operator {Hq} acts 
on V given by 

v=Y, ^(A)e* ek , 

kez L 

as 

({H Q }V)(A,0) = J2 ^ k ^k(AK 0k 

kez L 

A possible choice of T is 

(TV)(A,6)= 

We notice that this choice of T commutes with {-f^o}- 

The operator 1Z is the projector on the resonant part of the perturbation: 

nV = E ^(A)e' 0k = Vo(A), (5) 

wk=0 

since to is non-resonant. We also define the projector on the non-resonant part of the 
perturbation 

J\fV = ^2 V k {A)e iek = V - V . (6) 
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The global control follows directly from the definition of these operators T, 1Z and 
M : We construct a global control term for the perturbed Hamiltonian H + V, i.e. 
we construct / such that the controlled Hamiltonian H c = Ho + V + / is canonically 
conjugate to H + 1ZV. This conjugation is given by the following equation 

e^(H + V + f)=H + TZV, (7) 

where 

= E T^pW {rVT(nft + 1)V. (8) 

n=l 

We notice that if V is of order e, the control term / is of order e 2 . In general, the control 
term depends on all the variables A and 6, and acts globally on all phase space. 

Since oj is non-resonant, 1ZV = Vo(A) only depends on the actions and thus 
Ho + 1ZV is integrable. The derivation of Eqs. ©-© is given in Refs. j3E]. 

Starting from this global control, we derive a localised control such that the control 
term only acts in a given region of phase space around a selected invariant torus. We 
consider a nearly integrable Hamiltonian system : 

H(A,8) = H (A) + V(A,6). (9) 

We assume that H has the invariant torus with a non-resonant frequency vector u> at 
A = A . For V sufficiently small, the KAM theorem ensures that this invariant torus 
is preserved under suitable hypothesis. We expand Hamiltonian Q around A = A 
and we translate the actions such that the invariant torus with frequency uj is located 
at A = for Hq, and around A = for the perturbed Hamiltonian. Hamiltonian (|9j) 
becomes (up to a constant) 

H(A, 6) = uA + ev{0) + ew(0)A + Q(A, 6), (10) 

where Q is of order 0(||A|| 2 ), i.e., Q(O,0) = and d A Q(0, 0) = 0. Without any 
restriction, we assume that Hamiltonian (fTUj) is such that TZv = and TZw = : The 
mean value of v is absorbed into the total energy and the mean value of w into the 
frequency vector u>. 

For A sufficiently small, the perturbation 

V(A, 6) = ev(6) + ew(0)A + Q(A, 6), 

is small. We apply Eq. (jHJ) in order to get the control term /. However, for larger A, 
the control term is no longer small. Therefore we localise it in a region close to A = 0, 
i.e. we consider the following controlled Hamiltonian : 

H C (A, 0)=uA + V(A, 0) + /(A, ©Mil A||), 

where Q is a smooth characteristic function such that Q(x) = if x > 2e, and Q(x) = 1 
if x < e. The main drawback of this approach is that the control term is a priori of 
order e even if it is small since it is localised in a region near A = 0. In the next section 
we develop another approach where the control term / is of order e 2 and does no longer 
depend on A. 
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2.2. Localised control theory 

As in the previous section, we consider the family of Hamiltonians (|10jh For s = 0, 
Hamiltonian has an invariant torus with frequency vector uj located at A = 0. The 
problem of control we address is to slightly modify Hamiltonian (fTUj) near A = in the 
following way : 

H C (A, 0)=ujA + ev(0) + ew(0)A + Q(A, 0) + s 2 f(0)Q(\\A\\), (11) 

such that the invariant torus with frequency uj exists for the controlled Hamiltonian H c 
for higher values of the parameter e than in the uncontrolled case. Here Q denotes a 
smooth step function, meaning that the control only applies in a small part of phase 
space (of size s) : For instance, Q is a sufficiently smooth function such that Q(x) = 1 
if x < s and Q(x) = if x > 2s. Moreover, we notice that the control term / we apply 
is only a function of the angles in the region around the invariant torus. 

The main proposition of the localised control of Hamiltonian systems is the following 

one : 

Proposition 1: If v and w are sufficiently small and if v , w and Q are smooth, there 
exists a control term / such that the controlled Hamiltonian 

H C (A, 0)=ujA + sv(0) + sw(0)A + Q(A, 0) + s 2 f{0)tt(\\A\\), (12) 

is canonically conjugate to 

ff c (A,0) = £>A + Q(A,0), (13) 

where u> = uj + sa with a constant covector a and Q is of order 0(\\ A|| 2 ), i.e. Q(0, 6) = 
and <9a<3(0, 0) = 0. The control term / is given by 

/(0) = J\f [ w Tdv - s- 2 Q(-sTdv, 0)} , (14) 

which can also be written as 

f(0) = e - 2 M [V(0, 0)-V {-YdV{0, 0), 0)} . 

The important feature of the control term / is that it does only depend on the angle 
variables. Since the Hamiltonian H c has an invariant torus with frequency uj at A = 0, 
the controlled Hamiltonian H c has also this invariant torus in the region where A is 
close to 0. 

Proof: We consider the following transformations T mj b acting on functions V(A, 0) 

like 

(T m , b V)(A, 0) = e- md [V(e aA A + b, 0)] , 

where the covector m and the vector b are functions of from T L into 1R L . The operator 
m is the linear operator from R L to M. acting on a vector u G M. L as mu = mu. Here 
dm is the linear operator from R L to R L which is the product of the two linear operators 
acting on a vector u G 1R L as 

cMm = m'u + mu'. (15) 
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In Appendix A[ we check that the transformations T m b are canonical if b derives 
from a scalar function. We perform a transformation T m b on the controlled Hamiltonian 
H c given by Eq. (|12jl and we determine the functions m and b in the following ways : 

(i) The function b is determined such that the order e of the constant term in actions 
vanishes. 

(ii) The function m is determined such that the linear term in actions (which is of order 
e) vanishes. 

(Hi) The control term / is determined such that the constant term in actions [which is 
now of order e 2 after (i)] vanishes. 

We perform a transformation T em e b which is e-close to the identity : The expression 
of H c = T £mt£h H c is 

H c = e~ £md La + eub + ev(G) 

+ euJis! A + ewA + Q((l + en') A + eb, 0) 

+ e 2 w/i'A + e 2 wb + E 2 f(0)VL(\\(l + en') A + eb||)l , (16) 
where the covector n is defined by 



emd 



— 1 ^ (emd) n 



emd *-^(n + l)\ 

n=0 v ' 

The function \i is e-close to m : 

fi = m + smm' + ■ • • , 

and satisfies 

e e9ril A = A + A, 

so that 

1 + efi' = e £dlil ■ 1, (18) 

where e edm ■ 1 is a matrix, function of 9, which results from the action of the operator 
e e<9m on cons t an t function 1. First we notice that the scalar function 

q(A,0) = l-fi(||(l + sn')A + eb\\), 

is of order 0(||A|| 2 ). This can be seen from the equations 
g(O,0) = l-fi(e||b||), 

d A q(O,0) = -Q' ( e ||b||) +en'), 

where b is the transposed covector of b. The function h(0) will be chosen such that 
||b|| < 1 for all 0. Therefore we have q(0, 0) = and 8aq(O,0) = 0, according to the 
hypothesis on the smooth step function Q. 

Next, we expand the function Q((l + en') A + eb, 0) around A = : 



Q((l + en') A + eb, 0) = Q(eb, 0) + d A Q(eb, 0)(1 + en') A + Q 2 (A, 0), 
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where Q 2 is of order 0(||A|| 2 ), i.e. Q 2 (0,9) = and d A Q 2 (O,0) = 0. We notice that 
Q(eh,0) is of order e 2 and the covector d A Q(eb,6) is of order e since Q is of order 
0(||A|| 2 ). The Hamiltonian H c becomes 

H c = e~ £md La + eub + ev{0) 

+ eojfi' A + ewA + e 2 wfi' A + d A Q(eb, 0) (1 + £//) A 

+ Q(eb, 6) + e 2 wb + e 2 f(0) + Q 2 (A, 0) - f(d)q(A, 0)] . (19) 

The canonical transformation is determined by two equations : 

ub + v = 0, (20) 
Af [(w + £w + 9 A g( £ b, 6>))/*' + w + e- 1 9 A Q(eb, 6>)] = 0. (21) 

The control term is chosen such that 

f(0) = -Af [wb + e- 2 Q(eb, 6)] . 
Equations (J2DJ is solved in Fourier space. We expand the function b : 

b(0) = b k ei0k ' 

and the coefficients b^ are given according to Eq. (^0]) : 

, k 

b k = — T-v k , 

UK 

when ujk ^ 0, and bk = when uk = 0. Thus the vectorial function b is chosen to be 

b = -Tdv. 

We recall that we require that ||b|| = sup |b| < 1 which is ensured if v is sufficiently 
small and smooth and if u> satisfies a Diophantine condition (see the usual KAM proofs 
like for instance in Ref. [§]). In particular, we notice that this choice of b satisfies 
7Zb = 0. Equation (|2*Tj) is solved by choosing /x = Tp,, where p, satisfies 

fi + (ew + d A Q{-eVdv, 0))Tp.' + w + e^dxQi-eYdv, 6) = 0. (22) 

It is straightforward to check that Eq. (J2~T|) is satisfied if p, is a solution of Eq. ([22)1 . If 
the operator 1 + (ew + d A Q{— sYdv , 0))Td, which is e-close to the identity and then 
invertible, Eq. ()21j) has a solution 

/x = -r [1 + (ew + d A Q(-eTdv, 0))r9] _1 (w + e^dAQi-eVdv, 6)). 

The hypothesis on invertibility is fulfilled if w as well as v are small enough and smooth 
and if uj satisfies a Diophantine condition (see again Ref. |Hj). The covector /x which is 
e-close to m has the following expansion 



,x = _r w - e-Td A Q(-eTdv, 0) + 0{e). 



Localised control for non-resonant Hamiltonian systems 



9 



We notice that 7Z[i = since T and 1Z commute. The resulting Hamiltonian is then given 
by Eq. (JT3Jl where uj = u> + eTZfi, and Q = e~ em9 (Q2 — qf) which is of order 0(|| A|| 2 ). 
Note that we have dropped some additive constants of order e 2 , 1ZQ(—eTdv,6) and 
lZ{vfTdv ) since we recall that 1Z is the mean value with respect to the angles. The 
equations of motion for H c are 

6 = u + d A Q(A,0), 
A = -d e Q(A,6). 

Since Q(0, 6) = and <9aQ(0, 6) = 0, we see that A = is invariant, and that 
the evolution of the angles is linear in time with frequency vector u). Therefore, the 
Hamiltonian H c has an invariant torus located at A = with frequency vector uj. More 
precisely, the flow of the controlled Hamiltonian H c on A = is 

J{H c } T -l 
c I £m,eb 

where 

' e- e9ril (A - b) 
6 + em 






The equation of the torus A = is thus 




em - eb \ \ 6 + ein 



□ 



Remark 1 : We notice that if v — 0, the control term / given by Eq. (J 14)) is zero. 
In this case, the original Hamiltonian already has the invariant torus at A = 0. 

Remark 2 : Addition "property of the control term- In the case where more than 
one invariant torus needs to be created, we can add the control terms localised in non- 
overlapping regions of phase space. This is a straightforward extension to the previous 
case. The controlled Hamiltonian becomes 

M 

H(A, 6) = H (A) + eV(A, 0) + e 2 J2 M )^ (II A - A 2 ,||) , 

i=i 

where fi is defined for each region of phase space by Eq. (j!4|) . We notice that in each 
region of phase space, the operators T, 1Z and M are different since they are defined 
from the frequency vector of a given invariant torus. 
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3. Applications 

3.1. Forced pendulum Hamiltonian 

We consider the following forced pendulum model : 

H(p, x, t) = -p 2 + e [cos x + cos(x — t)] . (23) 

Figure H depicts a Poincare section of Hamiltonian (}23j) for e = 0.034. We notice that 
for e > 0.02759 there are no longer any invariant rotational (KAM) torus j7j. First, this 
Hamiltonian with 1.5 degrees of freedom is mapped into an autonomous Hamiltonian 
with two degrees of freedom by considering that t mod 2tc is an additional angle variable. 
We denote E its conjugate action. The autonomous Hamiltonian is 

H = E + y + s [cosx + cos(x - t)] . (24) 

The aim of the localised control is to modify locally Hamiltonian (|2*4*|) in order to 
reconstruct an invariant torus with frequency uj. We assume that u is sufficiently 
irrational in order to fulfill the hypotheses of the KAM theorem. First, the momentum 
p is shifted by uj in order to define a localised control in the region p w since the 
invariant torus is located near p w uj for Hamiltonian (}2*4*j) for e sufficiently small. The 
operators V and 1Z are defined from the integrable part of the Hamiltonian which is 
linear in the actions {E,p) : 

H Q (E,p) = E + oop, 



H = H Q + V, 



and Hamiltonian (j24|) is 
where 



p 2 

V(p, x, t) = e [cosx + cos(x — t)\ + —. (25) 
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The action of T, 1Z and TV on a function U of p, x and t given by 



are given by 



U(p,x,t) = U klM (p)e l(klX+k2t) ; 
(fci,fe 2 )ez 2 



ru= y Uk ^ e i(kix +k2 t) (26) 

(fcx,fc2)?t(0,0) V 1 ^ 

KU = U Q ,o(p), (27) 
A/"Z7 = ^ £4 llfe2 e i(feia;+fe2t) . (28) 

(fci,fc 2 )^(0,0) 



5.1.1. Global control The actions of T, 7?. and A/" on given by Eq. are 

iv = 



sin a; sin(x — t) 



on oj — 1 

A/V = £ [cos x + cos(x — t)] . 

Since TV depends only on x and t, and since V and are quadratic in p, it is 
straightforward to check that only the first two terms of the series (JHJ) are non-zero. The 
global control term reduces to 

/(p, a:, = -\{TV}{K + 1)V + ^{n/} 2 (2ft + l)V. 

Its explicit expression is given by 

fcosx cosfx — 1)\ e 2 ( cosx cosfa; — £)\ 2 , N 
/(p,a; > t) = ej) + 1 : - + TT + , • (29) 



a; u — 1 J 2 \ uj u — 1 

We notice that the control term is of order e, i.e. of the same order as the perturbation. 
However, the control term / acts only in a region where \p\ < e since it is multiplied 
by a function f2(|p|) such that fi(M) = 1 when \p\ < e, and ^(M) = when 
|p| > 2e. Consequently, the controlled Hamiltonian H (E, p) + V(p, x, t) + f(p, x, t)Q(p) 
is locally integrable (since it is locally conjugate to E+p 2 /2) provided that the canonical 
transformation is well defined (which is obtained when e is sufficiently small). A phase 
portrait of Hamiltonian (|23j) with the control term (^Sj) shows a very regular behaviour 
which persist for high values of e. However we notice that for e greater than one, the 
control term is no longer small compared with the perturbation. 

3.1.2. Localised control In order to apply the localised control as in Sec. 12.21 we notice 
that Hamiltonian (|2*%|) is of the form (fTT)J) with v = cosx + cos(x — t), w = and 
Q = p 2 /2. In this case the control term given by Eq. (JT3j) is equal to 

f^t) = - l -N(Td x v) 2 . 




Therefore the control term is equal to 

„. s 1 f cosx cos(x — 1)\ 2 1/1 1 \ , . 

This control term has four Fourier modes with frequency vectors (2, 0), (2, 2), (2, 1) and 
(0, 1). We consider the region in between the two primary resonances located around 
p = and p = 1. The control term given by Eq. ()30|) can be simplified by considering 
the region of phase space around p = 1/2. By keeping the main Fourier mode of this 
control term, i.e. the one with frequency vector (2, 1) which has the largest amplitude 
for lo close to 1/2, the control term becomes [3121 

fa(x,t) = 2 ^ (1 X _^ cos(2s-t). (31) 

For the numerical computations we have chosen lo = (3 — Vo)/2 (golden- mean invariant 
torus) which is the last invariant torus to break-up for Hamiltonian (|23|). 

A Poincare section of Hamiltonian ()23j) with the approximate control term (|3~T]) for 
£ = 0.06 shows that a lot of invariant tori are created with the addition of the control 
term precisely in the lower region of phase space where the localisation has been done 
(see Fig. EJ. Using the renormalization-group transformation we have checked the 
existence of the golden- mean invariant torus for the Hamiltonian H + e 2 f is given by 
Eq. (|3T)J) with e < 0.06965. By using the approximate and simpler control term f a given 
by Eq. (|3*T|) the existence of the invariant torus is obtained for e < 0.04857. However, we 
have checked using Laskar's frequency map analysis ^Hj that invariant tori and effective 
barriers to diffusion (broken tori) persist up to higher values of the parameter (e ps 0.2). 
The next step is to localize / given by Eq. (|3~U|) around a chosen invariant torus created 
by / : We assume that the controlled Hamiltonian H + e 2 f has an invariant torus with 
the frequency lo. We locate this invariant torus using frequency map analysis. Then we 
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Figure 3. (a) Poincare surface of section of Hamiltonian Ij23(l with the approximate 
control term (|32[1 with e = 0.034 (enlargement in the inset), (b) Same as (a) with 
initial conditions above the invariant torus. 



construct an approximation of the invariant torus of the Hamiltonian H + e 2 f of the 
form p = po(x,t). We consider the following localised control term : 

f^{p, x, t) = f(x, t)Q{\p - p (x, t)\), (32) 

where Q is a smooth function with finite support around zero. More precisely, we have 
chosen Q(x) = 1 for x < a, Q(x) = for x > (3 and a third order polynomial for 
x e]a, (3[ for which Q is a C 1 -function, i.e. Q(x) — 1 — (x — a) 2 (3/3 — a — 2x) /{(3 — a) 3 . 
The function p and the parameters a, (3 are determined numerically (a = 5 x 10~ 3 
and (3 = 1.5a). The support in momentum p of the localised control is of order 10~ 2 
compared with the support of the global control which is of order 1. 

Figure El shows that the phase space of the controlled Hamiltonian is very similar to 
the one of the uncontrolled Hamiltonian. We notice that there is in addition an isolated 
invariant torus. Using frequency map analysis [ID] , we check that this invariant torus 
corresponds to the one where the control term has been localised, i.e. its frequency is 
equal to (3 — Vo)/2. 

We notice that the perturbation has a norm (defined as the maximum of its 
amplitude) of 6.8 x 10~ 2 whereas the control term has a norm of 2.7 x 10~ 3 for e = 0.034. 
The control term is small (about 4% ) compared to the perturbation. We notice that 
there is also the possibility of reducing the amplitude of the control (by a factor larger 
than 2) and still get an invariant torus of the desired frequency for a perturbation 
parameter e significantly greater than the critical value in the absence of control. 

3.2. Delta-kicked rotor - standard map 

We consider the standard map : 

p n+ i = p n + K sin x n , (33) 
x n+ i = x n + p n+ i mod 27T. (34) 
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Figure 4. Phase portrait of the standard map for K = 1.2. 



This map is obtained by a Poincare section of the following Hamiltonian 

H(p, x, t) = — + — cos(x-mt). (35) 



Figure |U depicts a phase portrait of the standard map for K = 1.2. We notice that 
there are no KAM tori (dividing phase space) at this value of K (the critical value of 
the parameter for which all the KAM tori are broken is K c « 0.9716). Similarly to the 
forced pendulum, we consider the invariant torus with frequency u. By translating the 
momentum, we map Hamiltonian (J55j) into 



H = E + oop + e cos(x — mt) + — 
where e = Kj (47r 2 ). 



(36) 



m=—oc 



3.2.1. Global control Using the same computations as for the forced pendulum, the 
controlled Hamiltonian obtained by the procedure described in Sec. 12.21 is 

^2 



H(p, x, t) = ^— + e cos(x — mt) 



cos (a; — mt) e 



+oo 



uj — m 



cos(x — mt) 



2 \ z — ' u — m 

\m=—oo 



which is integrable and canonically conjugate to H = p 2 /2. We notice that although 
the control term is of the same order as the perturbation, it is more regular than the 
perturbation since its Fourier coefficients decrease like m _1 . 

3.2.2. Localised control Hamiltonian is of the form ([111)) with v = J2 m cos(x — mt) , 
w = and Q =p 2 /2. 
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Figure 5. Phase portrait of (a) the standard map for K = 5 and (6) the controlled 
standard map with the control term 



The control term given by Eq. (JT3j) becomes 

if cos(x — mt)\ 1 1 , 

/<*■«) = "2 ( E ^73^ + 4 E P7) 

\m=— oo / m=— oo v ' 

Again we notice that this control term is more regular than the perturbation since its 
Fourier coefficients decrease like m _1 . In particular, it is bounded in space and time, 
piecewise continuous in time. For instance, for uj — 1/2, the control term (}3T)) is equal 
to 

2 2 

f(x, t) = y cos(2s - t)x(t £ 2vrZ) + ^ X (t e 2vrZ). (38) 

The phase portrait of Hamiltonian (jHUj) with the control term (|3*H|) for K — 5 is 
depicted on Fig. We notice that in this case, the controlled kicked rotor is now 
a kicked pendulum : Instead of the rotor H = p 2 /2, the integrable part becomes a 
pendulum 

p 2 e 2 7T 2 
H = — + cos(2x - t), 

and the perturbation is a periodic 5-kick. We notice that the controlled Hamiltonian 
has invariant tori in the region near p = 1/2 (where the control has been localised). 
These invariant tori persist up to high values of the parameter K — ej (47r 2 ) larger than 
10 which has to be compared with K c w 0.97 in the absence of control. Note that the 
control term we use is bounded (conversely to the perturbation) and its amplitude is 
small compared with the amplitude of the Fourier coefficients of the perturbation (with 
a factor smaller than 10% depending on K). 

In order to recover a map, we need to locate the control at each u = (2m + l)/2 
for m G Z. For a given m G Z, the approximate control term is 

/i a) (M) = ^cos(2x-(2m + l)t). 
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Since each of these control term act on different regions of phase space, we sum all these 
control term to obtain a more global stabilization and recover a map. The control term 
becomes 



9 +00 
7T 



/( a )( x ,t) = _ ^ cos(2x - (2m + l)t). 



4 

m=— 00 

Next, we perform an inverse Fourier transform. The Hamiltonian becomes 

r 
I 



f( a \x,t) = — cos 2x [S(t - 2nm) - 5(t - (2m + . 



The controlled standard map is thus obtained by performing additional kicks: Each 
t = 2irm, in addition to the kicks of strength Kcosx, one has to perform kicks of 
strength (K 2 cos2x)/16, and each t = (2m + 1)tc, one has to perform kicks of strength 
— (K 2 cos2x)/16. This leads to the following form for the controlled standard map 

K 2 

P2m+1 = P2m + K Sin X 2m + — Sm2x 2m , 

lb 

K 2 . 

P2m+2 — P2m+1 ~ S m 2^2m+l) 

lb 



^2m+l — %2rn + 7:P2m+l? 



%2m+2 — %2m+l + -^P2m+2- 

It has a more compact form by using x n = x 2n and p n = p 2n '■ 

K 2 

Pn+i =p n + K sm x n + — sin 2x n 

16 

- sin(2x n + p n + Ksmx n + (K 2 sin2x n )/16)], (39) 
K 2 

x n+1 = x n +p n+1 + — sin(2x n +p n + Ksmx n + (K 2 sin2x n )/16). (40) 

If we neglect the kicks every t = (2m + 1)tc, we obtain the following map 

K 2 

Pn+i = Pn + Ksmx n + — sm2x n , (41) 

lb 

x n+1 = x n + p n+1 . (42) 

Conversely, if we neglect the kicks every t = 2mir, we get 

K 2 

Pn+i = Pn + K sm x n - — sm(2x n + p n + K sm x n ) , (43) 

Id 

K 2 

x n+ i = x n + p n+ i + — sm(2x n + p n + K sm x n ) . (44) 

A phase portrait of these maps are depicted on Figs. El Hand |H1 The most efficient 
control is obtained for the map (|39 )) -(|40 )) by considering the two additional kicks of order 
K 2 . The control term which only add the negative kicks does not lead to an efficient 
control although it appears slightly more regular than the uncontrolled case. Using 
frequency map analysis [TO], we have computed the critical thresholds for the break-up 
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X 



Figure 6. Phase portrait of the controlled standard map (|39(l - l|40(l for K = 1.2. 




A" 



Figure 7. Phase portrait of the controlled standard map l)41 [l -l)42 fl for K = 1.2. 

of the last KAM tori : There are invariant tori for the map (|4T)|) up to K c « 2.23 
which is more than twice the uncontrolled case (K c w 0.97). The map (|41| ) -(|42j l which 
is simpler than the map f|39j ) -(j40 p has invariant tori up to K c « 1.57. 

In this section, we have used the control for Hamiltonian flows in order to derive 
control terms for area-preserving maps. We note that a control method has been 
developed directly for area-preserving maps in Refs. fTTTj . 

3.3. Non-twist Hamiltonian 

We consider the following Hamiltonian 

H = E + up+ — + s (cos x + cos(x — t)) , (45) 
3 

where u = (a/5 — l)/2. A Poincare section of this Hamiltonian with e = 0.2 is depicted 
on Fig. The invariant torus with frequency u is located at p = for e = 0. We 
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Figure 8. Phase portrait of the controlled standard map (|43(l - (|44|l for K = 1.2. 



notice that this invariant torus is shearless since the second derivative of H with respect 
to p is zero on this torus. Hamiltonian (J4"5|) is of the form given by Eq. (jlO)) with 
v = cosx + cos(x — t), w = and Q = p 3 /3. The control term is given by Eq. (|14jl : 



The noticeable feature is that the modification of the Hamiltonian is of order e 3 
compared with the forced pendulum or the standard map where the control term is 
of order e 2 . A Poincare section of the controlled Hamiltonian (|4"5j) with the control term 
(j4fij) is depicted on Fig. El We notice that there are invariant tori in the region near 
p = that have been created with the addition of the small control term of order e 3 . 
For instance, for e = 0.2 and cu = (v5 — l)/2, the control term has a norm which is 
about 10% of the perturbation. In order to obtain a localised control, one has to apply 
the control term only in the region where the KAM invariant tori have been recreated, 
i.e., in the region p m 0. 
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Appendix A. The transformations T m b are canonical 

First we notice that T m b can be written as 



where T> is a dilatation and T a translation of the actions acting on a function V(A, 0) 

as 




(46) 



= e 




(D 1+m ,F)(A,0) = F((1 + a*')A,0), 



Localised control for non-resonant Hamiltonian systems 



19 



1.5 




X 



Figure 9. Poincare section for Hamiltonian 145JI with e = 0.2 and u = (V5 — l)/2. 




x 



Figure 10. Poincare section for Hamiltonian (|45|) with the control term l)46fl with 
e = 0.2 and u = (a/5 - l)/2. 



(W)(A,e) = y(A + b,0). 

A translation of the actions by a function b of 6 (i.e. independent of A) is obviously 
a canonical transformation if b derives from a scalar function, i.e. b = dj3. Thus it is 
sufficient to prove that T m o is also a canonical transformation. We notice that T m o is 
an automorphism since it is the product of two automorphisms : an exponential of a 
derivation and a dilatation. In what follows we prove that 

e- m9 D 1+M ' = e< mA >, (A.l) 

which is equivalent to say that T m o is a Lie transform generated by the scalar function 
mA. In other terms, Eq. (|A.1|) can also be written as [cf. Eq. JUJ)] 

2^™.! = e m9 e {mA} . 

Since these operators are automorphisms, it is sufficient to check that Eq. (|A.1J) is 
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satisfied on the basis (A, 6). First, we expand the operator e^ mA ^ : 



e {mA} _ e -ma+m'A0 A 



and we use the Trotter-Kato formula [T2j to express the exponential of the sum of two 
operators : 

e A+B = lim (e^e^)™. 



Any function of is invariant under the action of an operator e m ' A9 A/ n . Therefore it is 
straightforward to check that 



^ e -m5/n e m'Aa A /n^ n q = g-m9^ 



for all nGN. It follows that Eq. (jA.l|) is satisfied on 6. 
For A, we use the identity 

e MAd ANA = Ne M A ^ 



where M and TV are functions of 0. This identity follows from (M Ad^)NA = NMA. 
In particular, we have 

e A = l>„m . 



Using Eq. (|A.2|) we prove that 

^ e -ma/n e m'A0 A /n^ A _ e ~md / n ^m 1 / n A 

and it follows recursively that 

n 

^ e -m9/n e m'A9 A /n j ^ _ J~J g-fcmd/n ^m'/n ^ ^ 

k=l 

where the product is taken from right to left, i.e. YYk=i a k = a nd n -i ■ • ■ a>i- Concerning 
the operator e~ md Vi +fl i, the Trotter-Kato formula leads to the following expansion 



? X> 1+ , = lim e~ md (e m ' /n e md/n Y A, 

n—*oo \ J 



since 1 + (/,' — e 9m ■ 1 and drh = m' + md [see Eq. ()15|) ]. Using the same type of 
computations as for T m0 , we have 



o 

e m7n e m9/nV* _ TT ^kmd/n (^'/n 



n 

k=n— 1 

If we multiply the above expression by e~ md and change the index of the product 
{k' = n — k), it leads to 

e {mA} A = e - m9 £> 1+M ,A, 

and hence Eq. ()A.1|) is proved. 
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Remark : We notice that Eq. (jA.lj) implies that 

[L> e dih.i) — e u e -dzh.ie , 
which leads to the expression of the inverse of 1 + // : 

(e 9 ™ ■ l)' 1 = e m9 e- 9lh ■ 1. 
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